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Using extended mean field theory (EMFT) on the lattice, we study properties of the Higgs- 
Yukawa model as an approximation of the standard model Higgs sector, and the effect of higher 
dimension operators. We remark, as has been noted before, that the discussion of vacuum stability 
is completely modified in the presence of a (f>^ term, and that the Higgs mass no longer appears fine 
tuned. We also study the finite temperature transition. Without higher dimension operators the 
transition is found to be second order (crossover with gauge fields) for the experimental value of the 
Higgs mass Mj, = 125 GeV. By taking a (j)^ interaction in the Higgs potential as a proxy for a UV 
completion of the standard model, the transition becomes stronger and turns first order if the scale 
of new physics, i.e. the mass of the lightest mediator particle, is around 1.5 TeV. This implies that 
electroweak baryogenesis may be viable in models which introduce new particles around that scale. 


I. INTRODUCTION 

Even before the Higgs boson was discovered several stud¬ 
ies where the standard model (SM) couplings were run to 
very large energies were conducted. This resulted in both 
upper and lower bounds on the Higgs mass stemming 
from the triviality of the Higgs self-interaction and the 
electroweak (EW) vacuum stability, respectively. Sand¬ 
wiched between these bounds is a small region of Higgs 
masses for which the SM can be run at least up to the 
Planck scale mm- Due to the huge success of the SM 
when it comes to explaining results from accelerator ex¬ 
periments it was predicted that the Higgs mass would 
indeed lie inside this region and that the UV completion 
of the SM would enter only at the Planck scale where 
gravity becomes important. As the Atlas and CMS ex¬ 
periments announced the discovery of a Higgs-like parti¬ 
cle at 125 GeV, which is at least very close to the special 
region, these speculations about no new physics before 
the Planck scale got a lot of deserved attention. It is 
certainly interesting and important to thoroughly inves¬ 
tigate this possibility, see [1] and references therein. 

However, the apparent special value of the Higgs mass 
of course in no way excludes new physics at a lower scale 
and we will in this paper deal with generic aspects of UV 
completions at a fairly low scale of a few to tens of TeV. 
So let us depart from the arbitrarily postulated quartic 
self-interaction of the Higgs field. Since it does not in¬ 
clude dark matter or gravity, the SM is only an effective 
theory anyway, and there is no reason to assume a renor- 
malizable Higgs sector. We may add higher-dimension 
operators, which may not directly be relevant for the low 
energy physics but which can, for example, play an im¬ 
portant role for the stability of the EW vacuum or in the 
context of EW baryogenesis, in an effective field theory 
way. This will of course drastically change the running of 
the SM coupling constants so the question of whether the 
Higgs mass puts the universe in a near critical state or not 
loses its relevance in this context. We choose to consider 


generic higher-dimension operators instead of a specific 
UV extension of the SM because we can thus investigate 
aspects common to a broad class of models, which should 
make this study more relevant. It has previously been 
demonstrated [3H5], using the functional renormalization 
group (ERG) and various simplified versions of the SM, 
that the Higgs lower mass bound can indeed be lowered 
when higher dimension operators are included. 

Since the top Yukawa coupling is of order one, and be¬ 
cause we want to study the EW finite temperature tran¬ 
sition, it is desirable to use a nonperturbative approach, 
i.e. a lattice regularization of the model; for a perturba¬ 
tive study, see [5]. Unfortunately, it is not known how 
to regularize chiral fermion interactions on the lattice so 
we cannot consider the full SM. There are two sectors 
of the SM which can be studied separately, the gauge- 
Higgs sector, consisting of the weak gauge bosons and 
the Higgs field, and the Higgs-Yukawa sector, consisting 
of the Higgs field and the SM fermions. For a study of 
the first, see [THII]- Here, we choose the latter because 
of the large contribution to the SM Higgs sector from 
the Higgs-top interaction. The method we will employ 
to study this model is called extended mean field theory 
(EMFT), which has proven to be a highly accurate ap¬ 
proximation when applied to scalar field theories [HIS]. 
This allows us to obtain results with a computing effort 
orders of magnitude smaller than with full Monte Carlo 
simulations of the same model [HHsi, in addition to 
other advantages explained in the main body of this pa¬ 
per. See also m for a study of the Higgs-Yukawa model 
with an additional heavy fourth fermion family. 

The rest of the paper is organized as follows. In sec¬ 
tion |ll] we introduce the Higgs-Yukawa model and its lat¬ 
tice discretization followed by details on the implemen¬ 
tation of the chiral fermions in section mi Section m 
is devoted to the effective action in the EMFT approx¬ 
imation followed by our main results in section |V] We 
conclude in section ED 
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II. HIGGS-YUKAWA MODEL 


The Higgs-Yukawa model is a simplified version of the 
SM Higgs sector where the gauge degrees of freedom are 
neglected. The components of the model are the scalar 
complex doublet ip and the fermion doublets /. These 
couple to ip via Yukawa couplings with coupling constants 
Uf, which also determine the tree level fermions masses 
ruf = yfV after symmetry breaking, via the Higgs field 
expectation value (ip) = (0, Since the top quark 
is orders of magnitude heavier than the other fermions, 
it is common, and well justified, to restrict the fermion 
content to solely the top-bottom doublet dt = (t, by = 
itL,tiiyL,bRy ■ In these fields the Euclidean continuum 
action is given by 

[% d/, ^] = 5 h[(^] + [^, d/, , (1) 

with 

s'hM = : 

( 2 ) 

Sf[^, p] = jd'^x {^$d! + yb^Lpbn P yt^Lptn + h.c.} , 

( 3 ) 

where ip = iT 2 p\ T 2 is the second Pauli matrix and 
'^L/R = {tL/RibL/nY■ The other fermions (quarks and 
leptons) can, if desired, be added in a completely analo¬ 
gous way. 

In this study we will mainly be interested in the sym¬ 
metry broken phase and for notational convenience we 
will exploit the global SU(2) symmetry to make the ex¬ 
pectation value of p real and sit entirely in the lower 
component of p, i.e. we parametrize 

^ (°) ■ 

where v + h{x) is the Higgs field and gi{x) are the three 
Nambu-Goldstone modes. 

For the Higgs self-interaction we will consider higher 
dimension operators, in addition to the renormalizable 
interaction. The simplest extension is to add a dimension 
six contact term with six Higgs fields, {p''^py, as studied 
by ERG in [^, but operators of any dimension could just 
as well be included. Let us for now simply group these 
higher order terms in a “new physics” action m 


S'npM = 


d—6 i—1 


Ksm 


( 5 ) 


where is an operator of mass dimension d, rid is the 
number of operators with dimension d, Ci^d. are the Wil¬ 
son coefficients and Mbsm is the energy scale of the 
new physics, typically the mass of the lightest media¬ 
tor particle. Apart from p^ there is a second operator 


of dimension six, Oq = but it can be ne¬ 

glected if one assumes an approximate custodial sym¬ 
metry [19]. Naturally (naively), at low energy E, the 
effects of the higher-dimension operators are suppressed 
by factors {E/MbsmY~^- 


III. DIAGONALIZING THE OVERLAP 
OPERATOR FOR ARBITRARY CONSTANT 
HIGGS FIELD 


To efficiently, albeit approximately, integrate out the 
fermions we consider the Higgs field to be very slowly 
varying in space-time. Since we are primarily interested 
in the infrared properties of the model this assumption is 
reasonable. We thus assume that the fermions see a con¬ 
stant Higgs field, therefore the fermionic interaction can 
be diagonalized by going to Fourier space and the fermion 
determinant can be calculated without much effort. 

Due to the global SU(2) invariance, the fermion deter¬ 
minant can only depend on the magnitude of the Higgs 
field, \p^ = {v-\-hy+g\+g 2 -\-g‘^. Note that it depends on 
all the fields in p and not only on the expectation value. 
To simplify the derivation we apply an SU(2) transforma¬ 
tion such that p = (0, |i^|)^- Then, the different fermion 
flavors decouple and we have 

Sf = Y.N,j U^xjMff, ( 6 ) 

/ 

Mf = $+ yf\p\Ii, (7) 

where N^j is the number of colors for each fermion /, 
i.e. one for the leptons and three for the quarks. Unless 
otherwise specified we will include all SM fermions except 
the neutrinos with their Yukawa couplings set via the tree 
level relation yj = mfjv. 

For the model to be a realistic approximation to the 
SM Higgs sector it is important that the fermions be 
chiral. This is ensured by implementing the Neuberger 
overlap operator [20] when putting the fermions on the 
lattice. The overlap operator satisfies an exact lattice 
chiral symmetry which approaches the usual chiral sym¬ 
metry in the continuum limit a —>■ 0, a being the lattice 
spacing. Since we work with an effective model with a 
finite cutoff, this term will never completely go away but 
as long as the cutoff is well above the top mass the effects 
should be small. The overlap operator is given by 

jgiov) , A = 0 < p < 2r, 

a V yA^Aj a 

( 8 ) 

where is the usual Wilson operator with negative 

bare mass Mq and Wilson parameter r and p is a dimen¬ 
sionless parameter. In our calculations we will adopt the 
common choices ofr = I/2,p=I. The lattice action is 
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constructed by the following replacements: 
$ fL,RfR,L = fPR,LPR,Lf, 


( 9 ) 


Pr,l = 


I4±j5{h- 


= Pr.lT^D^°^\ 

2p 

( 10 ) 


after which the fermion operator, Eq. 0 , becomes 

^(OV) ^ ^(ov) ^ 1^1 . ( 11 ) 

In order to determine the fermion contribution to the 
action we will have to calculate the determinant, or 
equivalently the trace log, of this operator. This is most 
convenient in Fourier space where is diagonal. For a 
given 4-momentum p, its four eigenvalues come as com¬ 
plex conjugate pairs i'{p), (p), each with multiplicity 

two, where 


which can be calculated quite efficiently. Actually, since 
it only depends on one variable, y/ |(/7|, and will have to 
be evaluated very often, it will prove advantageous to 
precalculate it on a discrete set of values and interpolate 
to intermediate points. In summary, to a first approx¬ 
imation the effect of the fermions is the addition of an 
SU(2) symmetric contact term to the Higgs potential. 


It should be noted that the determinant is only posi¬ 
tive for generic Yukawa couplings in the approximation 
that the Higgs field is constant. Otherwise the Higgs 
field will fluctuate in the complex plane and introduce a 
sign problem, unless the fermions in each doublet have 
degenerate Yukawa couplings. 


r{p) = p 


L I i^/W+^f-p \ 

V + j 


f = = d^sin^ . 

/i ^ 


( 12 ) 

(13) 


Since the “mass term” pf \{p\ is real, the determinant of 

is real as well and the trace log takes the form of 
a real integral 


TrLog 



v{p) + pf \p\ 


(14) 

2p J ’ 


IV. THE EFFECTIVE ACTION AND EMFT 
SOLUTION 


For dehniteness we will consider only the \(j)\^ term in 
Snp, Eq. 0- Since we have no handle on the Wilson 
coefficient of this term we will set it to 1 and introduce 
Ae = (oMbsm)”^ where a is the lattice spacing. With the 
approximate treatment of the fermions above, we end up 
with the lattice action 



— K 


+ \px\^ + A {\ipxf' 


^WjTrLog 


'«l‘/5o|)) -kA6|^!>|^ 


(15) 


in terms of the conventional parameters as 

1 — 2A 

ap{x) = =-8 (16) 

At 

A = Xe = Sk^Xq. ^ _2^ ^ S^lS^p - AdKvho. (18) 

This action is quite similar to the complex model ^ 

studied by EMFT in [T3] and we can adopt the vector 
notation used there. 


= {hx,gi,x,g2,x,g3,xy + (D,0,0,0)'r = -f ($)'^, 

(17) 

with V = avjy/^, to derive the EMFT equations in an 
analogous fashion. Concentrating on the held at the ori¬ 
gin, 4>o, the hopping part of the action can be expressed 


The lattice without the origin is considered an external 
bath and will be self-consistently integrated out. This is 
equivalent to replacing the nearest-neighbor interaction 
term in Eq. (181 by its cumulant expansion with respect 
to the external bath. Truncating the expansion at second 
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order we obtain the following effective action: 

^EMFT = - A)$ + A (||$f - l)" 

+ ^iV,jTrLog(M(°") (y/V^||$|!)) (19) 

/ 

— 2v(v + h){2dK — Ai) + Ag ||^||^ , 

where A emulates propagation in the external bath. A 
closer inspection of the cumulant expansion reveals 

A = ^ 

±M.P 

= 2^2 diag((/i^/ip,5i,^gi,p,g2,A252,p,53,,x53,p))ext 

±M.P 

= diag(Ai, A 2 , A 2 , A 2 ), (20) 


where the diagonal form follows from the symmetries 
of the action. Since the action is still symmetric with 
respect to 0(3) rotations of gi,g 2 ,g 3 it is practical to 
rewrite the action in terms of two variables 


{v + ho) = 4>h, 



( 21 ) 

( 22 ) 


in which A = diag(Ai, A 2 ) and {gf^) = (<(>g) /3. The 
EMFT action is then given by 


^EMFT = (1 - Ai)<^^ + (1 - + X{^l + ^l-iy + Y iVc./TrLog (23) 

/ 

— 2v(j)h{2dK — Ai) + Ag {(fh + (/>g) , 


and the partition function becomes 


^EMFT =Af d(j)hd4>g 4>1 exp (—S'emft) : 


(24) 


where M is an irrelevant normalization constant. The 
unknown parameters v and A can then self-consistently 
be determined via the three self-consistency equations 


{M = (25) 

1 ^ y X- f V (26) 

(27r)'‘-k Ai - 2 kZ,, cos(pp) 

^_ (27) 

^ (27r)4 -k A 2 - 2KZh Ep cos(pp) ’ 

where {4>h)^ = {'Pfi) ~ {4’h)‘^ and the wave function renor¬ 
malization Zh is chosen such that the Nambu-Goldstone 
bosons are exactly massless m- The last two equations 
enforce that the connected 2-point function, from the ori¬ 
gin to the origin, is equal to its momentum-space ex¬ 
pression. The four- (or d-) dimensional integrals can be 
transformed into one-dimensional integrals by using the 
identity 


r d'^p 1 

J {2Tr)d a-Y)f, cos(pp) 


Jdre “^(/o(r))‘^, 
0 


(28) 


where Io{x) is a modified Bessel function of the first kind. 
See [13] for more details. 


A. Scale setting and observables 


In order to examine the effect of different cutoffs and 
possible higher dimension operators on the Higgs boson 
mass we need to set the scale of the lattice calculations. 
This is most naturally done by fixing the Higgs field ex¬ 
pectation value V to its phenomenological value. In terms 
of lattice variables we have 


y/^v 

V = - 

a 


= 246 GeV, 


(29) 


which, given v, determines the value of the cutoff A = 
1/a. Furthermore, we want to use physical quark masses, 
so we fix the Yukawa couplings using the tree level rela¬ 
tion iTif = yfV. We can now determine the Higgs boson 
mass Mfi in GeV as a function of the parameters of the 
Higgs potential by evaluating the inverse propagator at 
zero momentum: 




Ai - 2KZh Y cos(p^) 


KZh{{aMh)‘^ + (ap)"^), 

p—>-0 


( 30 ) 


M/ = 




-f Ai — SkZii 


_Y_ 

nZh 
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V. RESULTS 

In order to assess how well our EMFT method works in 
the presence of fermions, we compare it to already exist¬ 
ing full Monte Carlo results and results obtained using 
an analytic, approximate method, the constraint effective 
potential (CEP) [16]. In Eig. [^we show the Higgs expec¬ 
tation value in lattice units 0 as a function of the hopping 
parameter k at two different values of the new coupling, 
Ae = 0.1 in the upper panel and Xq = I.O in the lower 
one. Each color represents a different value of the quartic 
coupling A. In this comparison only the top and bottom 
quarks are included with degenerate Yukawa couplings 
Vb = yt = 175/246 and color factors = Nc,t = 1- 
It is clear that EMFT (solid lines) is a very good ap¬ 
proximation and gives results close to the Monte Carlo 
results (symbols) in all cases, in contrast to CEP (dot¬ 
ted lines) which works acceptably well for small values 
of Ae only. This is not surprising since the CEP calcu¬ 
lations in m are perturbative, whereas EMFT is fully 
nonperturbative. 

While we have no explanation for the remarkable ac¬ 
curacy of our EMFT approximation, we point out that 
a similar accuracy has been observed in the (p'^ lattice 
model near criticality (T^j: EMFT appears to perform 
well near a Gaussian critical point. 

For our actual results we will adopt a slightly differ¬ 
ent point of view on the model than the authors of m- 
Indeed, consider the origin of the higher dimension oper¬ 
ators. Generically they stem from UV completion of the 
SM and are thus associated with an energy scale which 
is typically the mass of the lightest of the “new” par¬ 
ticles which couple to the Higgs field. Let us call this 
energy scale Mbsm: where BSM stands for “beyond the 
standard model.” It is natural to assume that the cou¬ 
pling of this particle to the Higgs field is of order one, 
such that the coefficient in front of the \(j)f operator will 
be with the corresponding dimensionless coupling 

Ae = (oMbsm)”^- Notice, moreover, that A = a~^, not 
Mbsm, is the cutoff of the effective model, since it is di¬ 
rectly related to the maximum energy scale probed by 
the lattice action. In order to justify the effective treat¬ 
ment of particles heavier than Mbsm the cutoff A has to 
be sufficiently small. This leads to a hierarchy of scales 
condition, aMh I oMbsm, which in turn means 
Ae < I in the broken symmetry phase. As (oMbsm) is 
decreased toward I, more and more terms in the effective 
action would have to be taken into account in order to 
maintain a good approximation of the underlying model. 

It is quite challenging to preserve a good separation of 
scales while at the same time keeping the physical vol¬ 
ume large, and hence any lattice simulation is susceptible 
to large finite size effects. This is particularly true in a 
theory with massless modes, like the one we study here. 
In EMFT one generally works directly in the thermody¬ 
namic limit but in order to demonstrate the power law 
corrections coming from the Nambu-Goldstone modes, 
we have solved the self-consistency equations in a finite 



/t 



/t 


FIG. 1. The Higgs vacuum expectation value in lattice units 
as a function of the coupling k at various A, for Ae = 0.1 
(upper panel) and Ag = 1 (lower panel). The solid lines are 
EMFT calculations from this work and the pluses and dashed 
lines are full Monte Carlo simulations and CEP calculations 
respectively, both taken from [15]. The EMFT resnlts follow 
the Monte Carlo data closely for both values of Ag, whereas 
the CEP calculation gives reasonably accurate results for the 
upper, perturbative value only. 


volume. In Fig. we show the relative error on the Higgs 
mass as a function of the box size. The scale separation 
factor is aMbsm = a/IO (Ae = 0.1) for all values of a and 
we have marked both where the size of the correction 
is 50% of the mass itself and where an Ng = 32 lattice 
would be for two different lattice spacings. This demon¬ 
strates the need for very big lattices before one can even 
see the asymptotic power law scaling. 


A. The zero-temperature phase diagram 

It is most convenient to present the phase diagram in the 
(unphysical) bare parameters k, A and Ae. One can then 
pass to physical units via the Higgs expectation value in 
the broken phase. In Fig. [^one can see that the transi- 
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FIG. 2. Finite volume relative correction to the Higgs mass 
for two different lattice spacings calculated with EMFT with 
Ae = 1/10 and Mh = = 125 GeV. The horizontal line 

corresponds to a 50% correction. In order to see the asymp¬ 
totic corrections due to the massless Nambu-Goldstone 

modes, rather large lattices are needed, which poses a chal¬ 
lenge to full Monte Carlo simulations. In EMFT one can 
avoid the problem of thermodynamic extrapolation entirely 
by working directly in the thermodynamic limit. 

tion turns from second- to first-order as A is made more 
negative and in Fig. we show how the tricritical point 
depends on Ae. For A below the line the K-driven transi¬ 
tion is first order and above it is second order |2T]. Next, 
we fix Ag and look at the transition in the (A, /t)-plane. 
An example, where \q = 1/4, can be seen in Fig. 
The color of the line denotes the order of the transition, 
blue for second order and red for first order. The star 
marks the tricritical point and the arrow denotes how 
it moves as the number of lattice sites in the temporal 
direction Nt is decreased (see below for more details). 
The location of the tricritical point is obtained by fit¬ 
ting the critical vev on the first-order side with a power 
law (A) = c (Ac — A)^. There is a region close to the 
second-order line where one can perform calculations at 
a small lattice spacing a and since it is not possible to 
take the continuum limit of this effective theory anyway, 
it may also be viable to stay close to the transition on 
the first-order side. In fact, it turns out that a first or¬ 
der finite temperature transition will be found only there, 
denoted by the gray shaded area in the figure. 

B. Higgs mass lower bound 

Given a specific form of the Higgs potential, the lower 
bound on the Higgs mass is simply given by the mini¬ 
mal mass obtainable in the given parameter space of the 
potential [22]. For a pure potential it is found that 
the Higgs mass decreases when A decreases and the lower 
bound is thus obtained at vanishing quartic coupling. In 


FIG. 3. The tricritical point at zero temperature. For A below 
the line the transition is first order and above it is second 
order. 

0.118 
0.117 
0.116 
0.115 
0.114 
i; 0.113 
0.112 
0.111 
0.11 
0.109 
0.108 

-0.7 -0.75 -0.8 -0.85 -0.9 -0.95 

A 

FIG. 4. Zero- and finite-temperature transition in the (A, k)- 
plane at fixed Ae = 1/4. The transition turns from second 
to first order at the first-order endpoint marked by the star. 
The color of the line denotes the order of the transition, blue 
for second order and red for first order. As the lattice size is 
reduced in the temporal direction the endpoint moves along 
the arrow, and thus the gray shaded area marks the region 
in the plane where the finite temperature transition is first 
order. 


the Higgs-Yukawa model, this lower bound turns out to 
be just above 40 GeV and a positive A is needed to bring 
the Higgs mass up to 125 GeV. A negative coupling is 
obviously prohibited by the requirement of a bounded ac¬ 
tion. By introducing higher dimension operators, we can 
have a stable vacuum even at negative quartic coupling, 
and it is plausible that this could lead to an even lower 
Higgs mass. This was first demonstrated in |3], using the 
FRG on a chiral Z 2 Higgs-Yukawa model. Analogous re¬ 
sults were also obtained using the Higgs-Yukawa model 
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described above at physical values of the top and bottom 
masses in [4] and using the chiral Z 2 Higgs-Yukawa model 
plus an SU{3) gauge sector in [ 3 . Later, the authors 
of |16j came to similar conclusions using nonperturbative 
Monte Carlo simulations and perturbative CEP calcula¬ 
tions of the above described Higgs-Yukawa model with 
mass-degenerate top and bottom masses. Common to all 
these studies is that they add a (j)^ operator to the Higgs 
potential and when its coupling constant Xq is positive, 
the Higgs mass can be further reduced by making the 
quartic coupling A more and more negative. At some 
point, in the lattice regularized models, the phase transi¬ 
tion between the symmetric and broken phases turns first 
order (see Fig. and there is a hard lower bound on the 
lattice spacing a. Since the model is only effective this is 
in itself not a problem, but, since one wants aMh <C 1 , it 
bounds the region in parameter space where simulations 
are useful. 

Typically one finds that the Higgs mass goes to zero 
as one approaches the tricritical point from the second- 
order side although before zero is obtained one runs into 
subtle issues regarding new local minima of the effective 
action [53] . This means that the Higgs mass can in gen¬ 
eral be lowered by a large if not arbitrary amount, from 
its lower bound in the Ag = 0 case. This is demonstrated 
in Fig. where the Higgs mass is plotted as a func¬ 
tion of aMh for a few different values of uMbsm and A, 
together with the SM lower bound obtained at A = 0. 
As A is made more negative one approaches the regime 
of first-order transition and the Higgs mass decreases and 
can take values well below the SM lower bound. 



FIG. 5. The Higgs mass Mh as a function of aMh for various 
combinations of uMbsm and A, together with the SM lower 
bound on the Higgs mass obtained at Mbsm = cio and A = 0. 
In the presence of a operator a negative value of A is allowed 
and the Higgs mass can be lowered well below the SM bound, 
as indicated by the arrows. 

All in all it is clear that in the presence of higher di¬ 
mension operators the lower bound on the Higgs mass 
loses its meaning. In fact, one may argue that this is a 


null statement since the lower bound was calculated un¬ 
der the assumption that there is only the SM, and clearly 
new operators will change the running of the couplings. 


C. Finite temperature 

Let us now turn to the Mbsm dependence of the finite 
temperature transition. In gauge-Higgs systems it has 
been demonstrated that introducing a <p^ operator makes 
the phase transition stronger [siiiniiii], which in turn 
means that the critical mass, up to which the transi¬ 
tion is first order, increases. If it would increase past 
the observed Higgs boson mass, electroweak baryogene- 
sis might become possible again. Here, we present our 
findings using EMFT. Still in the infinite volume limit, 
we can introduce a nonzero temperature T = 1/ {aNt) by 
using a finite number Nt of lattice points in the tempo¬ 
ral direction. This gives us control over the temperature 
in discrete steps (for a fixed lattice spacing), so in order 
to get a good resolution one would need to work with 
rather fine lattices. However, this limits the range of 
available Mbsm because of the condition of scale sepa¬ 
ration oMbsm ^1- To overcome this problem we will 
determine the linear response of the system and then ex¬ 
trapolate to the desired temperature. Alternatively, to 
continuously vary the temperature, one could use a lat¬ 
tice action with anisotropic couplings. 

The observables of main interest are the critical Higgs 
mass for which the transition turns first order and the 
critical temperature. Another interesting observable is 
the Higgs mass for which the vev at the transition is of 
the same order as the critical temperature, which is the 
actual condition for a viable EW baryogenesis. All ob¬ 
servables depend on both Mbsm and the lattice spacing 
a so we need to calculate them in a two dimensional pa¬ 
rameter space. The starting point for determining this 
dependence is to obtain the phase diagram as in Fig. [^ 
for various oMbsm and Nt values, and then to determine 
how the first-order endpoint moves as a function of Nt- 
This is illustrated by the star and the arrow in Fig. [^ 
The finite temperature transition will be first order in 
the region between the T = 0 first-order line and the tra¬ 
jectory of the first-order endpoint, denoted by the shaded 
region in Fig. [^ We find that the endpoint moves on a 
straight line in the (k, A)-plane, as can be seen in Fig. 
for Ag = 1/4. By determining the lattice spacing a and 
the Higgs mass Mh at zero temperature along this line 
we can obtain the critical Higgs mass as a function of 
the BSM scale, shown in Fig. In Fig. we show how 
the critical Higgs mass for Ag = 1/4 changes for different 
strengths of the first order transition, measured in terms 
of (/c/Tc where t/c is the critical Higgs expectation value 
at the phase transition. In both of these figures the color 
of the line gives the critical temperature in GeV. 

Since our model is an effective one, we do not expect 
that the resulting curve is independent of the lattice spac¬ 
ing. Indeed, in the range of 2 < oMbsm < 4 we find 
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FIG. 6. The Nt dependence of the tricritical point for Ae = 
1/4. The trajectory is very well described by a straight line. 
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FIG. 7. The critical Higgs mass below which the finite tem¬ 
perature transition is first order for various values of oMbsm- 
The color coding gives the transition temperature in GeV. 
The 3 curves give a measure of the sensitivity of our effective 
theory to the cutoff. A higher cutoff makes the transition 
weaker. 


FIG. 8. The critical Higgs mass below which the finite tem¬ 
perature transition is first order for oMbsm = 2 and different 
strengths of the transition measured in terms of (j>clTc where 
(j>c is the critical Higgs vev at the transition. The color coding 
gives the transition temperature in GeV. 


This result is interesting since in many supersymmetric 
extensions of the SM one expects to find the lightest su¬ 
perpartners around this mass scale. The final scale sep¬ 
arations {aMh 1 ^ uATbsm) at the tricritical point 
for the different values of cMbsm are listed in Table |l] 
and it is evident that the validity of the effective model 
is somewhat strained, but not completely spoiled, due to 
the rather small separations. 


TABLE I. Separations of three scales in units of the inverse 
lattice spacing at the tricritical point for different values of 
olMbsm- 

Mbsm a~^ Mu 

2 1 0.19 

3 1 0.26 

4 1 0.29 


that the critical Higgs mass varies by about ten percent. 
This is in itself not a direct measure of the systematic 
uncertainties of the study since the choice of only includ¬ 
ing the operator introduces uncertainties which are 
hard to quantify, especially at the lower end of the inter¬ 
val. Moreover, different implementations of a UV cutoff 
will give somewhat different results. With these caveats, 
the variation of the critical line in the {Mh, MBSM)-plane 
gives a measure of the uncertainties within the model it¬ 
self. The final result is that for a Higgs mass of 125 GeV 
a BSM scale of around 1.5 TeV is needed to make the EW 
finite temperature phase transition first order, and this 
result changes only slightly even if we demand that the 
transition should be strongly first order with (pdTc ^ 1. 


Finally, to quantify the influence of our nonperturba- 
tive treatment of the fermions we repeated the calcula¬ 
tions considering just the Higgs sector and found that, 
after all, the fermions contribute only percent-level cor¬ 
rections to the purely bosonic case, see Fig. Remark¬ 
ably, the sign of the correction depends on the value of 
oMbsm- We also note that the BSM scale of « 1.5 TeV 
needed for a first-order finite-temperature transition is 
in good agreement with what one obtains in the gauge- 
Higgs model with a (/)®-term nnmi]. This could be used 
as an argument for leaving out the SM fermions and 
gauge fields from the simulations while studying higher 
dimension operators in the context of EW baryogenesis, 
and for including them perturbatively only. 
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Mbsm [GeV] 


FIG. 9. The critical Higgs mass below which the finite tem¬ 
perature transition is first order for various values of oMbsm 
with and without the SM fermions. The correction due to the 
fermions (shown by the arrows) is small and of indefinite sign. 

VI. CONCLUSIONS 

We have demonstrated that EMFT agrees very well with 
full Monte Carlo simulations of the Higgs-Yukawa model 
which contain both scalar fields and chiral fermions. We 
have then gone beyond what is possible for Monte Carlo 
simulations by lifting the mass degeneracy of the top and 
bottom quarks as well as including all other SM fermions. 
We have furthermore studied the EW finite temperature 


transition in the presence of a term in the Higgs poten¬ 
tial and thus obtained the critical Higgs mass for which 
this transition turns first order, something which has not 
been done nonperturbatively before. We find that with 
a BSM scale of about 1.5 TeV the transition turns first 
order for a Higgs mass of Mh = 125 GeV. At this point 
the effective model shows a separation of roughly a fac¬ 
tor 3 between the relevant scales Mh, a~^ and Mbsm- 
The value of Mbsm decreases only slightly if we demand 
a strong first order transition with (jic/Tc > 1. This 
scale is consistent with what is found in the gauge-Higgs 
model unmi], where the EW gauge fields are taken into 
account but the fermions are neglected, as well as in the 
perturbative study [6j. The scale is also only mildly de¬ 
pendent on the exact value of the lattice cutoff 1 /a within 
the window between Mh and Mbsm- It is however diffi¬ 
cult to assess the effect of neglecting other higher order 
operators. We have further shown that removing the 
fermions altogether shifts the critical Higgs mass only by 
a few percent, establishing that the Higgs sector itself 
is the dominating driving factor of the transition. To 
confirm that the gauge and fermion sectors always yield 
small contributions to the critical mass it would be inter¬ 
esting to study a model with different higher dimension 
operators, in particular one including both the Higgs field 
and the gauge fields. 
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